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We report on a study of jet shapes in inclusive jet production in pp collisions at y^s — 1.96 TeV 
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using the upgraded Collider Detector at Fermilab in Run II (CDF II) and based on an integrated 
luminosity of 170 pb _1 . Measurements are carried out on jets with rapidity 0.1 < |Y jct | < 0.7 and 
transverse momentum 37 GeV/c < P^ ot < 380 GeV/c. The jets have been corrected to the hadron 
level. The measured jet shapes are compared to leading-order QCD parton-shower Monte Carlo 
predictions as implemented in the PYTHIA and HERWIG programs. PYTHIA, tuned to describe 
the underlying event as measured in CDF Run I, provides a better description of the measured jet 
shapes than does PYTHIA or HERWIG with their default parameters. 

PACS numbers: 13.85.Ni, 13.85.Qk, 14.65.Ha, 87.18.Sn 



I. INTRODUCTION 



The measurement of the jet shape allows a study of the transition between a parton produced in a hard process and 
the collimatcd flow of hadrons observed experimentally yj. The internal structure of a jet is dominated by multi-gluon 
emissions from the primary outgoing parton and is expected to depend mainly on the type of parton, quark or gluon, 
creating the jet and the transverse momentum of the jet. In hadron-hadron collisions, the jet shape also receives 
contributions from initial-state radiation emitted from the colliding partons and multiple parton interactions between 
remnants (the so-called underlying event) . The effects of initial-state radiation are described by the parton showering 
in QCD Monte Carlo programs while the underlying event description is provided by phenomenological models. The 
comparison of jet cross section measurements with perturbative QCD predictions, as well as the estimation of QCD 
backgrounds in the search for new physics, requires an accurate description of the underlying event. The study of 
jet shapes at the Tevatron provides a precise means to test the validity of the models for parton cascades and the 
underlying event in hadron-hadron collisions. Measurements of the jet shape have been performed in pfj collisions at 

n n n 

i/s = 1.8 TeV |2 1, deeply inelastic scattering (DIS) 3| and photoproduction [j| processes in e ± p collisions at HERA, 
and e + e~ interactions at LEP1 5]. It was observed j| that the jets in pp collisions are significantly broader than 
those in e + e~ with most of the difference being explained in terms of the different mixtures of quark and gluon jets 
in the final state. The jets in DIS were found to be very similar to those in e + e~ interactions and narrower than 
those in pp collisions. In this paper, new jet shape results in pp collisions, based on CDF Run II data, are presented 
for central jets in a wide range of jet transverse momentum. For the first time, these measurements extend the study 
of jet internal structure to jets with transverse momentum up to 380 GeV/c. 



II. EXPERIMENTAL SETUP 



provides tracking over 



The CDF II detector is described in detail in 6J. In this section, the sub-detectors most relevant for this analysis 
are briefly discussed. As illustrated in Fig. ^ the detector has a charged particle tracking system immersed in a 
1.4 T magnetic field, aligned coaxially with the beam line. A silicon microstrip detector 
the radial range 1.35 to 28 cm. A 3.1 m long open-cell drift chamber, the Central Outer Tracker (COT) |8jL covers 
the radial range from 44 to 132 cm. The fiducial region of the silicon detector covers the pseudorapidity j^t range 
< 2, while the COT provides coverage for \r]\ < 1. The charged particles are reconstructed in the COT with 
a transverse-momentum resolution of <j(pt)/Pt ~ 1-7 1 10 _3 [GeV/c] _1 . Segmented sampling calorimeters, arranged 
in a projective tower geometry, surround the tracking system and measure the energy flow of interacting particles 
in \r)\ < 3.6. The CDF central barrel calorimeter lfj is unchanged from Run I and covers the region \rj\ < 1. It 
consists of an electromagnetic (CEM) calorimeter and an hadronic (CHA) calorimeter segmented into 480 towers of 

□ 



size 0.1 in rj and 15° in <fi. The end- wall hadronic (WHA) calorimeter 



complements the coverage of the central 



barrel calorimeter in the region 0.6 < \t]\ < 1.0 and provides additional forward coverage out to |ry| < 1.3. In Run 
II, new forward scintillator-plate calorimeters ^| replaced the original Run I gas calorimeter system. The new plug 
electromagnetic (PEM) calorimeter covers the region 1.1 < \rj\ < 3.6 while the new hadronic (PHA) calorimeter 
provides coverage in the 1.3 < \r)\ < 3.6 region. Each plug calorimeter is segmented into 480 towers with sizes that 
vary as a function of rj (0.1 in r\ and 7.5° in tft for \r]\ < 1.8 to 0.6 in r\ and 15° in (f> at I77I = 3.6). The calorimetry has 
a crack at r\ = (between the two halves of the central barrel calorimeter) and two cracks at rj = ±1.1 (in the region 
between the WHA and the plug calorimeters). The measured energy resolutions for electrons in the electromagnetic 
calorimeters are 14%/v^t (CEM) and 16%/%/E© 1% (PEM) where the units are expressed in GeV. The single-pion 
energy resolutions in the hadronic calorimeters, as determined with test-beam data, are 75%/\/Et (CHA), 80%/y/E 
(WHA) and 80%/Ve<S)5% (PHA). Cherenkov counters located in the 3.7 < \rj\ < 4.7 region [13j measure the average 



number of inelastic pp collisions per bunch crossing and thereby determine the beam luminosity. Finally, a three-level 
trigger system |l4| is used to select events online, as described in section V. 



III. MONTE CARLO SIMULATION 



Monte Carlo event samples are used to determine the response of the detector and the correction factors to the 

n 

hadron level [15J for the measured jet shapes. The generated samples are passed through a full CDF detector 



n 



simulation (based on GEANT3 [lfj where the GFLASH |l7j package is used to simulate the energy deposition in the 

e analysis chain as in the dai 
3 and HERWIG 6.4 Q 



calorimeters) , and then reconstructed and analyzed using the same ^analysis chain as in the data. Samples of simulated 
inclusive jet events have been generated using the PYTHIA 6.203 |18j and HERWIG 6.4 [19| Monte Carlo generators. 
In both programs, the partonic interactions are generated using leading-order QCD matrix elements, including initial- 
and final-state parton showers. CTEQ5L ^cj parton distribution functions are used for the proton and antiproton. 
The HERWIG samples have been generated using default parameters. The PYTHIA samples have been created 



using a special tuned set of parameters, denoted as PYTHIA- Tune A 



21 1, that includes enhanced contributions from 



initial-state gluon radiation and secondary parton interactions between remnants. Tune A was determined as a result 
of dedicated studies of the underlying event in dijet events performed using the CDF Run I data 2^. In addition, 
two different PYTHIA samples have been generated using the default parameters with and without the contribution 
from multiple parton interactions (MPI) between the proton and antiproton remnants. The latter are denoted as 
PYTHIA- (no MPI). The HERWIG samples do not include multiple parton interactions. Fragmentation into hadrons 



is carried out using the string model 
model IJ in HERWIG. 



23 as implemented in JETSET 24] in the case of PYTHIA and the cluster 



IV. JET RECONSTRUCTION 



An iterative cone-based midpoint algorithm |26j in the Y-(f> plane |9j is used to reconstruct jets from the energy 
deposits in the calorimeter towers for both data and the Monte Carlo simulated events, and from final-state particles 
for the Monte Carlo generated events. This procedure is explained in detail below for the jet reconstruction from 
the calorimeter towers. In the first step, the electromagnetic and hadronic sections of each calorimeter tower are 
preclustered into a physics tower. The position of each section is determined from the unit vector joining the vertex 
of the interaction and the section's geometrical center. Each section is assumed to be massless. The four-vector 
components of each physics tower are then computed using the four-momentum sum of its electromagnetic and 
hadronic sections; only towers with transverse momentum above 0.1 GeV/c are further considered. In a second step, 
each physics tower with transverse momentum above 1 GeV/c is used to define a seed for the jet search. Starting 
from the seed with highest transverse momentum, a cone is drawn around each seed and the physics towers inside a 



distance V( Ar ) 2 + (A0) 2 < R/2, with R = 0.7 are used to determine the direction of the new cluster as indicated 



9 



in Eqs.[DandH 



£ C luster = J2 E towc \ P^ st ™ = ^ i = X,y,Z (1) 



phys. towers phys. towers 



i pcluster i pcluster pcluster 
pcluste, = A / (p c lu8ter)2 + (p c lust e r)2; F c luS te r = ___±^ i _ ) fluster = ^-1 ( ^__ ) (2) 



where y cluster and ^ clustcr denote the rapidity and azimuthal angle of the cluster, respectively. Starting from the list 
of resulting clusters, the procedure is iterated until the contents of the clusters remain unchanged. In a third step, the 
midpoint (Y-(f> plane) between each pair of stable clusters separated by less than 2R is added to the list of clusters. 
The clustering algorithm, as explained above, is again iterated until stability is achieved. This latter step gives the 
name to the jet algorithm and was introduced in order to address the theoretical difficulties |2Jj of the cone-based 
jet algorithm used in Run I 28| . Finally, the cone size is expanded from R/2 to R \^(\ and the momentum sharing 
of overlapping clusters is considered. Overlapping jets are merged if their shared momentum is larger than 75% of 
the jet with smaller transverse momentum; otherwise two jets are formed and the common towers are assigned to 
the nearest jet. The variables for jets reconstructed from the calorimeter towers are denoted by f-pcAL' ^cal anc ^ 
4>cal- As mentioned, the same jet algorithm is applied to the final-state hadrons in Monte Carlo generated events. 
In this case, the four-vector components of each individual hadron are used as input to the algorithm and no cut on 
the minimum transverse momentum of the particles is applied. The variables of the hadron-level jets are denoted by 

pjot ^jot , ,jet 

T.HAD > 'HAD duu ™HAD 1 

The reconstruction of the jet variables in the calorimeter is studied using Monte Carlo event samples and matched 
pair of jets at the calorimeter and hadron levels. These studies indicate that the angular variables of the jet, Iqal 
and <^cal' are reconstructed in the calorimeter with no significant systematic shift and with a resolution, for jets 
with P^cal > 20 GeV/c, of the order of 0.02 units and 0.025 units, respectively. The resolutions improve as the 
measured jet transverse momentum increases. The jet transverse momentum measured in the calorimeter, P^cali 
systematically underestimates that of the hadron level jet. This is mainly due to the non-compensating nature of the 
calorimeter 29) . For jets with -BJ?q AL > 20 GeV/c the jet transverse momentum is reconstructed with an average 
shift of —20% and an r.m.s of 17%. The reconstruction of the jet transverse momentum improves as -Rj^cal increases. 
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For jets with P^cal > 130 GeV/c the jet transverse momentum is reconstructed with an average shift of —12% and 
an r.m.s of 9%. An average correction is extracted from the Monte Carlo using the following procedure: matched 
pairs of jets are used to study the difference between the jet transverse momentum at the hadron level, f-THAD' an( l 
the corresponding measurement in the calorimeter, -P-^cal' The resulting correlation is used to extract multiplicative 
correction factors, C(P^ t CAh ), which are then applied to the measured jets to obtain the corrected jet transverse 
momenta, Pj^cor. — C x -Ptcal l^oj - 



V. EVENT SELECTION 



This analysis is based on a sample of inclusive jet events selected from the CDF Run II data corresponding to a 
total integrated luminosity of 170 pb _1 . Events were collected online using three-level trigger paths, based on the 
measured energy deposits in the calorimeter towers, with several different thresholds on the jet transverse energies. In 
the first-level trigger, a single trigger tower with transverse energy above 5 GeV or 10 GeV, depending on the trigger 
path, is required. In the second-level trigger, a hardware-based clustering is carried out where calorimeter clusters are 
formed around the selected trigger towers. The events are required to have at least one second-level trigger cluster 
with transverse energy above a given threshold, which varies between 15 and 90 GeV for the different trigger paths. 
In the third-level trigger, jets are reconstructed using the CDF Run I cone algorithm |28| and the events are required 
to have at least one jet with transverse energy above 20 to 100 GeV depending on the trigger path. Offline, jets are 
reconstructed using the midpoint algorithm, as explained above, starting from seed calorimeter towers with transverse 
momentum above 1 GeV/c and only considering towers with a minimum transverse momentum of 100 MeV/c in the 
clustering procedure. The following selection criteria have been imposed: 

• One reconstructed primary vertex with z-component, Vz, in the region \Vz\ < 60 cm. Events with more than 
one primary vertex are removed to eliminate contributions from pile-up events with multiple proton-antiproton 
interactions per beam crossing. The data used in this study was collected at Tevatron instantaneous luminosities 
in the range between 0.2 x 10 31 cm _2 s _1 and 4 x 10 31 cm _2 s _1 for which, on average, less than one interaction 
per crossing is expected. 

• tf T /VEr < 3.5 GeV 1/2 , where tf T (E T ) denotes the missing (total) transverse energy of the event as determined 
from the energy deposits in the calorimeter towers. This cut eliminates beam-related backgrounds, beam halo 
and beam-gas contributions, and cosmic rays. 
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• At least one jet with P^ COK > 37 GeV/c and Y^ ct in the region 0.1 < |y jot | < 0.7. 

The cut on the minimum P^cor ^ s dictated by the trigger. In order to avoid any possible bias on the measured jet 
shapes due to the three-level trigger selection, the thresholds on P^cor' applied to the different data samples, have 
been selected such that the trigger is fully efficient in the whole kinematic region under study. The measurements are 
performed for central jets in a rapidity region away from calorimeter cracks and inside the fiducial region of the CDF 
tracking system. 

VI. JET SHAPE 
A. Jet shape definition 

The differential jet shape as a function of the distance r = \J AY 2 + A(f> 2 to the jet axis, p(r), is defined as the 
average fraction of the jet transverse momentum that lies inside an annulus of inner radius r — Sr/2 and outer radius 
r + Sr/2 around the jet: 



jets 

where 7Vj et denotes the total number of jets, Pr(r — Sr/2,r + 5r/2) is the transverse momentum within an annulus 
and the jet shape is determined for values of r between 0.05 and 0.65 using 5r = 0.1 intervals. The points from the 
differential jet shape at different r values are correlated since, by definition, p(r) 5r = 1. 

The integrated jet shape, *(r), is defined as the average fraction of the jet transverse momentum that lies inside a 
cone of radius r concentric to the jet cone: 



tt (r) = J_\-^4, 0<r<R (4) 

where, by definition, \&(r = R) = 1. The integrated jet shape is determined in intervals 5r = 0.1 between r = and 
r = 0.7, and the points at different r values are strongly correlated. 
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B. Jet shape reconstruction 

Calorimeter towers are used for both data and Monte Carlo simulated events to reconstruct the differential jet shape. 
For each jet, the scalar sum of the transverse momentum of the calorimeter towers assigned to it, PT(r — Sr/2, r+Sr/2), 
with a distance to the jet axis r' = y/(Y towor - yJ ct ) 2 + (0t°wcr — 0J ct ) 2 between r - Sr/2 and r + Sr/2, is determined 
and divided by Pt(0,R)- The differential jet shape, p CAL (r), is then determined following the prescription in Eq. 3. 
Similarly, the integrated jet shape, ^f CAh (r), is reconstructed using the calorimeter towers as defined in Eq. 4. The 
same procedure is applied to the final-state particles in Monte Carlo generated events to reconstruct the differential 
and integrated jet shapes at the hadron level, Py^f(r) an< ^ ^mcPM' respectively. In the case of hadron-level jets no 
grid in the (Y-<fi) space has been used. 



C. Jet shape using charged particles 

The CDF tracking system provides an alternative method to measure the shape of the jets using charged particles. 
For each jet, tracks with transverse momentum, p^ ack , above 0.5 GeV/c and pseudorapidity, rf 1 ' , in the region 
j^trackj ^ 4 are assigned to it if their distances, r, with respect to the jet axis are smaller than 0.7, and the tracks 
project to within 2 cm of the z-position of the primary vertex. The differential and integrated jet shapes, p TRK (r) 
and *™ K S(r), are then reconstructed using the track information and following Eqs. 3 and 4. The measured jet 
shapes using tracks are employed to study systematic uncertainties on the central measurements as determined 
using calorimeter towers (see next section). Therefore, detailed studies have been performed on track reconstruction 
efficiency inside jets as a function of r and the jet and track transverse momenta, for both data and simulated events, 
using track embedding techniques [^|. The difference between efficiencies in the data and Monte Carlo are about 3%, 
and approximately independent of r for tracks with 0.5 GeV/c < p^ ack < 2.0 GeV/c. For tracks with p^ ack > 2.0 
GeV/c, the difference in efficiency is of the order of 5% at the core of the jet, decreasing as r increases up to r = 0.5. 
For r > 0.5 no difference in efficiency is observed. The effect on the reconstructed jet shapes is smaller than 0.5% and 
thus has been absorbed into the systematic uncertainty. 



VII. UNFOLDING AND SYSTEMATIC STUDIES 

The measured jet shapes, as determined using calorimeter towers, are corrected back to the hadron level using 
Monte Carlo samples of generated events. PYTHIA-Tune A provides a good description of the measured jet shapes 
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in all regions of P^ and is used to determine the correction factors in the unfolding procedure. 



A. Jet shape corrections 

The measured jet shapes are corrected for acceptance and smearing effects back to the hadron level. The correction 
factors also account for the efficiency of the selection criteria and for jet reconstruction in the calorimeter. Differential 
and integrated jet shapes are reconstructed with Monte Carlo samples using both calorimeter towers, p^C ( r ) an d 
f^ L (r), and final-state particles, PMC D ( r ) ano - ^mcPMi m different regions of Py COH and -R|?had' respectively. 
Correction factors, defined as D(r) — Pmc° ( r ) I Puc ( r ) ana - ^( r ) = ^mc D ( r ) /^MC ( r ) > are then com P u ted separately 
in each bin of P^cor- The corrected differential and integrated measurements are determined from the measured jet 
shapes as p(r) = D{r) ■ p CAL (r) and 9(r) = I{r) ■ \& CAL (r). The correction factors D(r) do not show a significant 
dependence on Pjf and vary between 1.2 and 0.9 as r increases. For the integrated jet shapes, the correction factors 
7(r) differ from unity by less than 10% for r > 0.2. 



B. Systematic uncertainties 

A detailed study of the different sources of systematic uncertainties on the measured jet shapes has been per- 
formed 



3(]: 



• The measured jet transverse momentum has been varied by ±5% in the data to account for the uncertainty on 
the determination of the absolute energy scale in the calorimeter. The effect on the measured jet shapes is of 
the order of 2%. 

• The unfolding procedure has been repeated using bin-by-bin correction factors extracted from HERWIG instead 
of PYTHIA-Tune A to account for any possible dependence on the modeling of parton cascades. The effect on 
the measured jet shapes is about 2% to 5%. 

• The ratios of uncorrected jet shape measurements as determined using calorimeter towers and tracks, 
p (r)/p T (r) and ^ rCAL (r)/^ rTRKS (r), are compared between data and Monte Carlo simulated events. 
The deviations from unity observed in the data/Monte Carlo double ratio, below 5% for the whole P^ et range, 
are included in the systematic uncertainties to account for the uncertainty on the description of the inactive 
material in front of the calorimeter and its response to low-energy particles. 
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• The measurements are performed in different periods of Tevatron instantaneous luminosity (between 0.2 • 
10 31 cm _2 s _1 and 4 • 10 31 cm _2 s _1 ) to account for possible remaining contributions from pile-up events. No 
significant effect is found. 

The total systematic uncertainties on p(r) and ^(r) have been computed for the different r ranges by adding in 
quadrature the deviations from the central values. The statistical uncertainties are negligible compared to the sys- 
tematic errors except for jets with P^ ct > 300 GeV/c. The systematic uncertainties have been added in quadrature to 
the statistical errors and the total uncertainties are shown in the figures. The total uncertainty in the measured data 
points, for different P^ ct and r ranges, varies between 5% to 10% except for jets with Pj? * > 300 GeV/c for which the 
total error is above 20%. 

VIII. RESULTS 

The corrected differential and integrated jet shapes, p(r) and \&(r), refer to midpoint jets at the hadron level with 
cone size R = 0.7 in the region 0.1 < |Y jct | < 0.7 and 37 GeV/c < Pjf* < 380 GeV/c. 

A. Comparison with Monte Carlo 

Figures[21and|21show the measured differential jet shapes, p(r/R), in bins of P^ ct for jets in the region 0.1 < |Y> et | < 
0.7 and 37 GeV/c < Pf ot < 380 GeV/c, compared to the PYTHIA-Tune A and HERWIG Monte Carlo predictions 
at the hadron level. The measured jet shapes show a prominent peak at low r which indicates that the majority of 
the jet momentum is concentrated at r/R < 0.2. At low P^ et , the fraction of transverse momentum at the core of the 
jet is about a factor of 6 times larger than that at the tail. This factor increases at higher Pj, et and is of the order 
of 100 for jets with Pj?* > 340 GeV/c. PYTHIA-Tune A provides a good description of the measured jet shapes 
in all regions of Pj?\ The jets predicted by HERWIG follow the measurements but tend to be narrower than the 
data at low P^ et . The latter can be attributed to the absence of additional soft contributions from multiple parton 
interactions in HERWIG, which are particularly important at low Py Ct . 

Figures E| and present the measured integrated jet shapes, ^S(r/R), in bins of Pjf', for jets with 0.1 < |yj ot | < 0.7 
and 37 GeV/c < P^ < 380 GeV/c, compared to HERWIG, PYTHIA-Tune A, PYTHIA and PYTHIA-(no MPI) 
predictions, to illustrate the importance of a proper modeling of soft-gluon radiation in describing the measured jet 
shapes. 
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Figure shows, for a fixed radius tq = 0.3, the average fraction of the jet transverse momentum outside r — tq, 
(1— ^(tq/R)), as a function of Pjf . The points are located at the weighted mean in each Pj?* range. The measurements 
show that the fraction of jet transverse momentum inside a given fixed r^/R increases (1 — ^(tq/R) decreases) with 
P^ et , indicating that the jets become narrower as P^ t increases. PYTHIA with default parameters produces jets 
systematically narrower than the data in the whole region in PJ?'. The contribution from secondary parton interactions 
between remnants to the predicted jet shapes (as shown by the difference between PYTHIA and PYTHIA- (no MPI) 
predictions) is important at low Pj?\ PYTHIA- Tune A predictions describe all of the data well (a % 2 test in Fig. El 
gives a value of 13.6 for a total of 18 data points). HERWIG describes the measured jet shapes well but produces jets 
slightly narrower than the data at low Pjf . This results in a significantly higher \ 2 value of 33.8 for 18 data points. 

B. Quark- and gluon-jet contributions 

Figures [7| and |H1 present the measured integrated jet shapes, ^(r/R), in bins of Pf ct , for jets with 0.1 < \Y^ ct \ < 0.7 
and 37 GeV/c < P^ ot < 380 GeV/c, compared to PYTHIA- Tune A predictions (as in Figs. 0]andEJ). In these figures, 
predictions are also shown separately for quark- and gluon-jets. Each hadron-level jet from PYTHIA is classified as 
a quark- or gluon-jet by matching (Y-<f> plane) its direction with that of one of the outgoing partons from the hard 
interaction. The Monte Carlo predictions indicate that, for the jets used in this analysis, the measured jet shapes 
are dominated by contributions from gluon-initiated jets at low Pjf* while contributions from quark-initiated jets 
become important at high P^ et . This can be explained in terms of the different partonic contents in the proton and 
antiproton contributing to the low- and high-Py Ct regions, since the mixture of gluon- and quark-jet in the final state 
partially reflects the nature of the incoming partons that participate in the hard interaction. Figure [!|] shows the 
measured 1 — ^(tq/R), = 0.3, as a function of P^ ct compared to PYTHIA- Tune A predictions with quark- and 
gluon-jets shown separately. The trend with P^ et in the measured jet shapes is mainly attributed to the different 
quark- and gluon-jet mixture in the final state and perturbative QCD effects related to the running of the strong 
coupling, a s (P^ et ) |5j. The Monte Carlo predicts that the fraction of gluon-initiated jets decreases from about 73 % 
at low P^ et to 20 % at very high P^ et , while the fraction of quark-initiated jets increases. 
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IX. SUMMARY AND CONCLUSIONS 

Jet shapes have been measured in inclusive jet production in pp collisions for jets in the kinematic region 37 GeV/c < 
i-ijf* < 380 GeV/c and 0.1 < |Y jct | < 0.7. Jets become narrower as P^ et increases which can be mainly attributed to the 
change in the quark- and gluon-jet mixture in the final state and the running of the strong coupling with P^ et . PYTHIA 
Monte Carlo predictions, using default parameters, do not give a good description of the measured jet shapes in the 
entire * range. PYTHIA- Tunc A, which includes enhanced contributions from initial-state gluon radiation and 
secondary parton interactions between remnants, describes the data better. HERWIG gives a reasonable description 
of the measured jet shapes but tends to produce jets that are too narrow at low P^ et which can be attributed to the 
absence of soft contributions from multiple parton interactions in HERWIG. Jet shape measurements thus can be 
used to introduce strong constraints on phenomenological models describing soft-gluon radiation and the underlying 
event in hadron-hadron interactions. 
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FIG. 1: Longitudinal view of half of the CDF II detector. 
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FIG. 2: The measured differential jet shape, p(r/R), in inclusive jet production for jets with 0.1 < |y Jct | < 0.7 and 37 GeV/c < 
Pj? < 148 GeV/c, is shown in different regions. Error bars indicate the statistical and systematic uncertainties added in 
quadrature. The predictions of PYTHIA-Tune A (solid lines) and HERWIG (dashed lines) are shown for comparison. 
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FIG. 3: The measured difFerential jet shape, p(r/R), in inclusive jet production for jets with 0.1 < |Y Jct | < 0.7 and 148 GeV/c < 
< 380 GeV/c, is shown in different P^f regions. Error bars indicate the statistical and systematic uncertainties added in 
quadrature. The predictions of PYTHIA-Tune A (solid lines) and HERWIG (dashed lines) are shown for comparison. 
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FIG. 4: The measured integrated jet shape, ^(r/R), in inclusive jet production for jets with 0.1 < |Y Jct | < 0.7 and 37 GeV/c < 
Py Ct < 148 GeV/c, is shown in different Py Ct regions. Error bars indicate the statistical and systematic uncertainties added 
in quadrature. The predictions of PYTHIA- Tune A (solid lines), PYTHIA (dashed-dotted lines), PYTHIA-(no MPI) (dotted 
lines) and HERWIG (dashed lines) are shown for comparison. 
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FIG. 5: The measured integrated jet shape, ^(r/i?), in inclusive jet production for jets with 0.1 < |F Jct j < 0.7 and 148 GeV/c < 
Py Ct < 380 GeV/c, is shown in different regions. Error bars indicate the statistical and systematic uncertainties added 
in quadrature. The predictions of PYTHIA- Tune A (solid lines), PYTHIA (dashed-dotted lines), PYTHIA-(no MPI) (dotted 
lines) and HERWIG (dashed lines) are shown for comparison. 
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FIG. 6: The measured 1 - *(0.3/i?) as a function of P^ for jets with 0.1 < |Y jct < 0.7 and 37 GeV/c < P^ < 380 GeV/c. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. The predictions of PYTHIA- Tune A (solid 
line) , PYTHIA (dashed-dotted line), PYTHIA-(no MPI) (dotted line) and HERWIG (dashed line) are shown for comparison. 



25 



m data 0A < | y*| < 0.7 

- PYTHIA Tune A 



1 



^ 0.5 



0.25 


' 1 



7^0.75 
^ 0.5 



01 



0.25 


1 1 



>0.75 
^ 0.5 
0.25 








•* nh (on — Ipt 

— - - UIUUI 1 JC l 

2_ .-y/' ■■■■ quark-jet 

: / P T M > 37 GeV/c 
P T * < 45 GeV/c 

i i i i 1 i i i i 1 


P T * < 55 GeV/c 

1 1 1 1 1 I 1 1 I 1 


P T iet < 63 GeV/c 

i i i i 1 i i i i 1 


P T iet < 73 GeV/c 

I T I I 1 1 I I T 1 


*'* 

— ** f* 

— / * 

// 

_ f 

I ' P T M > 73 GeV/c 
P T * < 84 GeV/c 

1 T I I 1 I T T I 1 


■" 

•' / / 

/ 1 
/ / 
#/ 

~~ / 

P T * > 84 GeV/c 
P T M < 97 GeV/c 

I I T I 1 T I I T 1 


— •* 
.* 

.* 

** /' 

— ' / * 

— ft 

— J / 

— •/ 

— / 

: P T Jet > 97 GeV/c 
P* < 1 1 2 GeV/c 

i i i i 1 i i i i 1 


* y"' 

— .* // 

— •* / / 

•/ 

/ 

" P* > 1 1 2 GeV/c 
P T * < 1 28 GeV/c 

i i i i 1 i i i i 1 


• f ' 

' / * 
* / * 

~ Mi 

— 1 
/ 

: P T iet > 1 28 GeV/c 
P T Jet < 1 48 GeV/c 

i i i i 1 i i i i 1 



0.5 



1 



0.5 



1 



0.5 



r/R 



r/R 



r/R 



FIG. 7: The measured integrated jet shape, ^(r/R), in inclusive jet production for jets with 0.1 < Y Jct < 0.7 and 37 GeV/c < 
Py Ct < 148 GeV/c, is shown in different regions. Error bars indicate the statistical and systematic uncertainties added 
in quadrature. The predictions of PYTHIA- Tune A (solid lines) and the separate predictions for quark-initiated jets (dotted 
lines) and gluon-initiated jets (dashed lines) are shown for comparison. 
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FIG. 8: The measured integrated jet shape, ^(r/i?), in inclusive jet production for jets with 0.1 < j F Jct j < 0.7 and 148 GeV/c < 
Py Ct < 380 GeV/c, is shown in different regions. Error bars indicate the statistical and systematic uncertainties added 
in quadrature. The predictions of PYTHIA- Tune A (solid lines) and the separate predictions for quark-initiated jets (dotted 
lines) and gluon-initiated jets (dashed lines) are shown for comparison. 
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FIG. 9: The measured 1 - 9(0.3/ R) as a function of P^ for jets with 0.1 < |Y jct | < 0.7 and 37 GeV/c < < 380 GeV/c. 
Error bars indicate the statistical and systematic uncertainties added in quadrature. The predictions of PYTHIA- Tune A 
(solid line) and the separate predictions for quark-initiated jets (dotted line) and gluon-initiated jets (dashed line) are shown 
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